The magnetic states of lightly hole-doped cuprates in the clean limit. 
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We have performed extensive zero field /^SR experiments on pure YBa2Cu306+j, and diluted Y- 
rare-earth substituted Yo^Euo.osI^CuaOe+j/ and Yo.925Ndo.o75Ba2Cu306+ H at light hole-doping. 
A common magnetic behavior is detected for all the three families, demonstrating negligible effects 
of the isovalent Y-substituent disorder. Two distinct regimes are identified, separated by a crossover, 
whose origin is attributed to the concurrent thermal activation of spin and charge degrees of freedom: 
a thermally activated and a re-entrant antiferromagnetic regime. The peculiar temperature and 
hole density dependence of the magnetic moment m(h, T) fit a model with a (spin) activation 
energy for the crossover between the two regimes throughout the entire investigated range. The 
magnetic moment is suppressed by a simple dilution mechanism both in the re-entrant regime 
(0 < h < 0.056) and in the so-called Cluster Spin Glass state coexisting with superconductivity 
(0.056 < h < 0.08). We argue a common magnetic ground state for these two doping regions and dub 
it frozen antiferromagnet. Conversely either frustration or finite-size effects prevail in the thermally 
activated antiferromagnetic state, that vanishes at the same concentration where superconductivity 
emerges, suggesting the presence of a quantum critical point at h c = 0.056(2). 



I. INTRODUCTION 

The magnetic ground state of cuprates at very low dop- 
ing is qualitatively understood as that of a frustrated 1 
two dimensional Heisenberg antiferromagnet (2DHAF). 
The 2DHAF nature was recognized very early 2 - from the 
successful description that its long wavelength approxi- 
mation, the nonlinear a model (NLcrM), provides of the 
finite correlation lengths detected by neutron scattering 3 
above the transition to three dimensional (3D) order. 

The wave function of a hole in the Neel AF background 
leads to a local spin singled, i.e. it neutralizes a Cu spin, 
but its dynamical effects reduce the ordering tempera- 
ture T/v rather more effectively than a mere magnetic- 
site dilution. The abrupt reduction of the ordering tem- 
perature T^v (as detected very early on by implanted 
muons 5 ) corresponds also to a drop of the static mag- 
netic moment^— at Q = (hh)- The strong frustrating 
effect of holes in the 2DHAF background may be under- 
stood in several ways. Staggered spin spirals^ have been 
predicted and their hole density dependence reproduces 
the incommensurate magnetic scatteringiS detected in 
La2- a; Sr a ;Cu04. Stripes have also been predicted, either 
as charged spin solitonsii or from electronically driven 
phase separation^ They have been detected^ 3 - in spe- 
cific doping and distortive conditions (e.g. La1.6-xNdo.4- 
Sr^CuC^), and are possibly present also as correlated 
fluctuations. 

Notwithstanding this partially quantitative agreement 
between models and experiments we are still far from a 
complete understanding of the low-doping Mott Hubbard 
insulator and the number of open questions grows on 
approaching and then crossing the metal-insulator (MI) 
transition. They concern the precise nature of the mag- 
netic ground state versus temperature, the origin of the 
spin glass phase intervening at slightly higher doping, 



the presence of static magnetic moments in the very low- 
doping superconducting state, and ultimately the com- 
petitive or cooperative role of magnetic interactions in 
the strongly correlated superconductor. 

The theoretical difficulties increase with hole density, 
but the scarcity of experimental studies in the less puz- 
zling insulating range prevents precise checks across the 
whole insulating phase diagram. In addition the real 
cuprate compounds, as opposed to the idealized CuC>2 
model layer, are already distinguished at very low doping 
by structure and intrinsic disorder.— It is now clear,— ~— 
although not always recognized, that YE^CuaOe+j/, 
where the charge reservoir is farther removed from the 
Cu02 layer, represents the closest physical system to 
the clean limit, whereas randomly placed divalent cations 
make La2- £! ;Sr a ;Cu04 and Yi_. r Ca, / Ba2Cu306 dirty limit 
cases, where disorder, increasing with doping, actively 
modifies the scenario. 

In this work we decided to investigate extensively the 
low-doping range, starting from the clean limit case. 
Besides pure YBa2Cu306+ y (Y100% in the following) 
we addressed two additional compositions, to investi- 
gate the effect of small lattice distortion on the magnetic 
phase diagram by substitution of Y with a rare earth: 
Y .92Euo.o8Ba 2 Cu30 6+ j / (Eu8%), yielding a very small 
mismatch of the cation radii, and Yo.92sNdo.o75Ba2Cu3- 
Oq+ v (Nd7.5%), yielding equivalent mismatch to that of 
the 5% Ca substitution investigated in a parallel work: 20 
It turns out that isovalent substitution affects negligibly 
both the magnetic behavior and the MI transition. 

Muon spin spectroscopy (/j,SR) was selected as a lo- 
cal probe sensitive to static magnetic order, even with 
small moments, and to local spin rearrangements. This 
technique reveals a number of evident trends, yielding 
quantitative constraints for a model of doped 2DHAF 
appropriate for clean cuprates. We carefully identify two 
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different magnetic regimes and measure their distinctive 
parameters, evidencing their criticai behavior with dop- 
ing in the whole investigated range. 

The paper is organized in three main sections: Sec. |TT] 
briefly reviews models and experimental facts; Sec. IIIII 
describes the muon results in terms of a model of ther- 
mally activated charge excitations; Sec. IIVI discusses the 
implication of the data and of the model on the mag- 
netic state in doped Mott-Hubbard insulating cuprates. 
The appendices describe the relation between the mag- 
netic field at the muon site and the average staggered Cu 
moment, they provide details on the samples and on the 
muon data analysis. 



II. THE DOPED S = \ 2D HEISENBERG 
ANTIFERROMAGNET. 

Copper perovskites, characterized by one or more 
Cu02 layers of square plaquettes, have a common start- 
ing point in the two dimensional (2D) Neel order of their 
parent compound, technically a half-filled Mott-Hubbard 
insulator of the charge transfer— type, dominated by a 
large isotropic antiferromagnetic (AF) exchange J < 
along the in-plane bonds, with a much reduced coupling 
J' along the c axis. 

The predicted transition for a cubic system with the 
same J along all the three directions (3D) would be 
rather large, Xjvo ~ 1500 K, whereas the observed tran- 
sitions are T N = 310,420 K for La 2 Cu0 4 and YBa 2 - 
CU3O6, respectively, and very large in-plane correlation 
lengths 6 persist far above Tjy. This behavior, previously 
encountered in localized spin systems, such as the fluorite 
perovskites^ is described at high temperature by the 2D 
Heisenberg (2DHAF) model. For a pure 2D system the 
Mermin- Wagner theorem 2 ^ predicts Tjy = 0, but three- 
dimensional order sets in at an intermediate temperature, 
scaling roughly as l/log(J/J'). The ideal 2DHAF be- 
havior with its large quantum spin reduction for S = 1/2 
and the high temperature correlations are quantitatively 
described 2 in the NLerM long-wavelength approximation. 
The spin stiffness p s of cuprates places them in the renor- 
malized classical regime, where, even for J' = 0, a Neel 
state would be present at T — 0, as opposed to the quan- 
tum disordered regime, where fluctuations would kill the 
order parameter even at zero temperature. 

The Heisenberg exchange interaction is obtained from 
the Hubbard model at half-filling in the t-J approach. 
Departing from half-filling, i.e. doping the parent com- 
pound with a holes density h, the measured magnetic 
transition Tjsf(h) drops rapidly until a low temperature 
inhomogeneous magnetic order is established, often re- 
ferred to as cluster spin glass (CSG), with onset temper- 
ature Tf. The CSG phase is the close precursor of the 
metal insulator transition, taking place at a critical den- 
sity h s where superconductivity arises. It is a precursor 
both in the trivial sense that it lies in between the AF and 
the SC phases, and because increasing the temperature 



at any small finite density h leads to metallic transport 
properties^ The qualitative understanding of the mag- 
netic hole-doped ground state is that these holes localize 
on oxygen^ at low temperatures and develop a large ef- 
fective ferromagnetic (F) exchange K > 0. The strong 
reduction of the transition, dTjy/dh, is then the effect of 
a magnetic frustration^ A hole in the AF background ac- 
tually produces a current around the plaquette, that can 
be mapped to a static spin twist 2 ^ and still described, 
in the hydrodynamic limit, by the NLctM, as it is done, 
e.g., in the spiral models^^ Alternatively it may lead to 
diagonal or parallel stripes*^ 

The glassy nature of the CSG stated is however 
controversial, since magnetic structure coexists with 
disorde r 10 ' 27 and the observed feature s 28 ' 29 are not ex- 
clusive of a true spin glass state. We simply acknowledge 
that increasing hole density, hence frustration, leads to 
a more disordered state, whose experimental signature 
is the abrupt flattening of the magnetic transition de- 
pendence on hole density, i.e. dTf/dh <C dT^/dh. In 
all cuprates this so-called CSG phase survives well inside 
the superconducting dome, in the form of low tempera- 
ture static magnetis m 14 ' 16 ' 30 " — embedded in the metallic 
phase. The coexistence is intrinsic since it takes place 
over a range of hole densities h — h s < 0.02 — 0.03, typi- 
cally ten times greater than the inevitable hole inhomo- 
geneity of stoichiometric origin. 

A peculiar feature of the AF low doping region is 
the crossover between two thermal regimes, identified 
by the re-entrant behavior of the magnetic moment 
in La2- a; Sr a ;Cu04, YBa2Cu306+j/ and Yi-^Cayl^Cua- 
06- 139 La Nuclear Quadrupole Resonance (NQR) 31 
and muon o 14 ' 16 ' 31 ' 32 ' 34 ' 35 identify a thermally activated 
regime where the static moment is strongly reduced pro- 
portionally to hole density and a re-entrant regime where 
the moment recovers a nearly hole-independent zero tem- 
perature value. The crossover temperatures qualitatively 
agree with the onset of spin dynamics freezing, detected 
e.g. by Cu NQR 3 - 6 - and muonS ) 14 ' 34 but also with the acti- 
vation temperatures of the variable range hopping regime 
in transport j 15 ' 37 ' 38 Neutrons observe a decreasing static 
moment at (| h I) in the re-entrant state, associated with 
an increase of the rod (5 \ q) elastic intensity 6 (I integer, 
q real) that matches the muon moment recovery. It is 
worth stressing that in these conditions it is improper to 
consider only the moment deduced from the Bragg peak 
intensity at (^ as if a simple Neel state were present. 
Partial and full sum rulesS 9 - are hard to implement and 
imply very large errorbars, whereas local probes offer a 
more direct, almost model independent way of determin- 
ing the full static moment. 

Besides these commonalities, remarkable differences 
are experimentally detected among the cuprates. We 
focus here on a few relevant points, limiting ourselves 
to the cuprate families stemming from the two parents 
La2Cu04 and YBa2Cu30g (in this order, unless other- 
wise specified, for the values in parentheses): i) Distinct 
Tjy; ii) Distinct critical densities for the suppression of 
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FIG. 1: (color on-line) Temperature dependence for three Y- 
Ba2Cu3C>6+i/ samples of a) the magnetization (Eq. [TJ and 
b) the magnetic volume fraction defined in Sec. [C] Vertical 
arrows indicate the magnetic transition. 

T N (h c « 0.02,0.06); Hi) A region of pure CSG behav- 
ior, without superconductivity, sizeable^ in La2_ x Sr x - 
CuC>4, almost vanishing in YBa2Cu3C>6+j,; iv) Distinct 
critical densities for the onset of the superconducting T c 
(h s sa 0.05,0.06). The latter are rather close to each 
other, but clearly different^ as it is confirmed also by 
Ca doping in YBa 2 Cu306+ y , leading to h s = 0.071 for 
8% Ca, see Ref. Hl|. These main points originate from 
the strong influence of the disordered heterovalent sub- 
stituents (Sr 2+ , Ca 2+ ), which place L^-zSr^CuC^ at 
the dirty limit and YBa2Cu306+ y at clean limiti 17 ' 18 

III. EXPERIMENTAL RESULTS. 

Magnetic moments are proportional to the spin S and 
mostly reside on Cu ions in the square plaquettes. In 
YBa2Cu306+y, as discussed in details in Appendix \A\ 
the internal field at the apical oxygen (AO) muon site 
provides a measure of the staggered moment m(h, T) ac- 
cording to 

m(h, T) = (i2nB AO (h, T)/B AO (0), (1) 

where fj,2D ~ 0.6 Bohr magnetons, is the value of the 
magnetic moment appropriate to the spin 1/2 2DHAF, 
with its quantum spin reduction* 2 - and B A o(0) is the 
muon internal field value at zero temperature in the un- 
doped compound. This determination is independent 
of the details of the magnetic structure, provided that 
it preserves a local in-plane, two-sublattice collinear ar- 
rangement (Appendix [5]) . Conservation of this local spin 
order condition allows muons to detect also the short 
range CSG state. 



Figure QJi, shows the behavior of m for three YBa2- 
CU3O6+J; samples whose hole density ranges from nearly 
undoped, to just over the critical value where supercon- 
ductivity appears, the strongly underdoped regime. Fig- 
ure^ displays the fraction of implanted muons detecting 
the magnetically ordered state, extracted from the muon 
asymmetry analysis (details in Appendix [C]) , which mea- 
sures the fraction of the sample volume belonging to the 
magnetic state. 

First of all, the data of the two panels allow the deter- 
mination of the magnetic transition temperature T m (we 
do not distinguish yet between Tjy and Ty, because the 
muon cannot directly identify the two states). Secondly, 
at intermediate doping the order parameter m deviates 
from a standard power law behavior like that of the un- 
doped sample (dashed lines). The power law, with lower 
T m and rescaled magnetization, is followed only at high 
temperature, whereas an upturn (solid line) appears to- 
wards the undoped zero temperature value, m(0) = 1. 

For finite h the recovery of the undoped magnetization 
m(0) (Fig. []J,) justifies the term re-entrance for this be- 
havior. A similar upturn is just barely perceptible also 
in the nearly undoped case of Fig. [T^,; it has been de- 
scribed before, first^^i in La 2 _ K Sr a; Cu04, then 3 ^ 3 ^ 
in YBa2Cu306+ y and many more instances are shown 
in Fig. [U Muon and NQR agree quantitatively on the 
low temperature upturn, granting that this is not one of 
those very rare instances where the muon alters its local 
surroundings. 

The data are naturally described as the crossover from 
a low temperature magnetization towards a magnetically 
weaker high temperature regime, whose origin must be 
found in the thermal activation of charge and spin de- 
grees of freedom. The functional form of the upturn 
looks like an exponential decay to the rescaled behav- 
ior. Thus, in order to distinguish the two regimes, we 
perform a preliminary fit of the data to the function 34 
m(h,T) = [m A + (m F - m A ) e - cT ] (1 - T/T m f, where 
fnR.A stand for the moment in the re-entrant and acti- 
vated regimes, respectively. Notice that the former corre- 
sponds to the experimental ratio at T — 0, mp — m(h, 0) 
and the latter, m A , yields the zero temperature extrap- 
olation indicated by the dashed line in Fig. [TJi. For 
h > 0.055 we must revert to a simple power law fit 
m(h,T) = mp (1 — T/T m )^. Similar fits are obtained 
for all samples, including those with Eu and Nd partial 
substitutions. 

A compact view of the transition temperatures, with 
the moment as a color map is provided by Fig. [2] for 
all these fits, for the three families of cuprates under 
investigation. It is evident that they display a unique 
behavior, within experimental errors (by the way, this 
observation is remarkably confirmed in all the data pre- 
sented in this section). Another straightforward consid- 
eration is that T m agrees qualitatively with many pre- 
vious observations^ ^ 16 ' 43 on unsubstituted YBa2Cu3- 
Oe+y, but only here it is presented as a function of cali- 
brated hole density h. We recognize two regions: the first 




FIG. 2: (color on-line) Common phase diagram of YBa2Cu3- 
06+y, Y .92Euo.o8Ba2Cu306+ H and Yo. 9 2sNdo.o75Ba2Cu3- 
Oe+y' transition temperatures Tm and Tf from the fits to 
Eq. [2] experimental moment m from fit and Eq. [1] encoded 
in the nonlinear color mapping (see right bar); hole density 
h, see Sec. [B] for details. The phase boundaries are best fits 
with h c = h s = 0.056(2) (see text) and Ta labels the guide 
for the eye of Fig. Inset: zoom of the superconducting T c , 
from magnetic susceptibility. 



for < h < 0.05 where T m drops rapidly with increasing 
hole density, and we identify it with the Neel transition; 
the second for h > 0.05, where the T m (h) dependence is 
much weaker, and we associate this transition with Tf. 
The Neel temperature follows a quadratic behavior vs. 
hole density T N0 (1 - (h/h c ) 2 ) with T N0 = 422(5) K and 
critical concentration h c — 0.056(2) as shown by the solid 
line in Fig. [2] 

The color map used in Fig. [2] is nonlinear to empha- 
size the upturn of m(T) (Fig. [I}, identifying a crossover 
between the two regimes. The re-entrant region is charac- 
terized by a uniform dark red color towards T = 0, signi- 
fying that rap — m[h, 0) approaches the full 2D moment 
H2D at all dopings, even underneath the superconducting 
dome. This is shown by the parameters m,F,A plotted in 
Fig. [3^ (crosses). The re-entrant zero temperature pa- 
rameter mp is very close to one, up to the highest hole 
densities, until it undergoes an abrupt reduction around a 
distinct h c p s» 0.08. Conversely the activated parameter 
tua, plotted in the same figure, displays a clear distinct 
critical behavior, fitted to a power law tua = (1 — h/h c ) a 
with a = 0.35(5) and h c = 0.058(2), equal within error- 
bars to the value of Fig. [3] 

It is tempting to include the wiAih) power- law behav- 
ior of Fig. directly into the fit model, identifying its 
doping dependence as m^(/i) = mp(h) (1 — h/h c ) a . It 
is however clear that a constant hole density does not 
reproduce the thermally activated features of m and a 
straightforward extension is to replace h by Ha{T) — 
h exp(— Ta/T), assuming for the sake of the argument 
that delocalized (hopping) holes are responsible of the 




FIG. 3: (color on-line) Dependence on hole doping of: a) Re- 
entrant moment tuf, and activated moment tua from the fit 
to Eq.[2](phis and cross, same quantities from the preliminary 
fit, see text); b) Width of the moment distribution measured 
by muons above and below Ta- Squares for Y100%, circles 
for Eu8% and triangles for Nd7.5%. 

magnetization reduction. This leads to 
m(h,T) = m A (h,T) (l - 

This activated-hole function yields excellent best fits 
of the data as it is shown by the solid lines in Fig. [4] 
We note however that Eq. [5] could also hold in a stripe 
scenarioii, where it is possible to imagine a mechanism 
by which the influence of thermally activated excitations 
on magnetism depends on the stripe density, which in 
turns is proportional 44 to h. The difference between Eq.[2] 
and the functional form of the preliminary fit is often 
beyond resolution. The critical temperatures (T/v and 
Tf) fitted in the two models are almost identical, so that 
Fig. [2] is unchanged. The new parameter Ta, is shown 
in Fig. [5] and its guide for the eye is reported also in 
Fig. [2] as a quantification of the crossover temperature 
(dashed line). The activation energy follows a dome that 
replicates roughly the AF phase boundary, albeit with a 
large down-shift. On entering the CSG phase coexisting 
with superconductivity, for h > h c , the fit can still be 
performed, but Ta is hardly distinguishable^ from Tf, 
hence we report only the latter in Fig. [21 

The two moments displayed in Fig. [3^, have a coun- 
terpart in the activated- hole model (Eq. [2]): the zero 
temperature moment, mp, is still a free parameter, 
whose functional dependence on h is directly deter- 
mined by best fit procedures, whereas tua is obtained as 
limr-^oo m,A(h, T). These new determinations (colored 
symbols) agree perfectly with the previous ones (crosses) 
as it is shown in the figure. 
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FIG. 4: (color on-line) Temperature dependence of the mo- 
ment m with the best fit to Eq. [2] for the three families. 



Two very distinct critical behaviors are also shown in 
panel b, where we plot the widths Attif.a of the two 
static moment distributions. These widths are obtained 
directly from the static relaxation of the muon asymme- 
try precession as <7ao{T)/2-k"{Bao{T) (Sec. [C]), taking 
cuts through the data at two temperatures, the lowest 
one, 2 K, below all Ta and another one, 70 K, above the 
largest Ta- Both quantities display a divergence which 
points, within errors, to the same two distinct critical 
densities mentioned above, h c and h C F- 




FIG. 5: (color on-line) Dependence on hole doping of the 
activation temperature Ta from the best fits to Eq. [2] The 
line is a guide for the eye and the data for h > he represent 
the transition temperature Tf. 



the hole concentration h = 1/8 and by the shift of the 
MI transition in heavily Y-substituted- 7 - YBa2Cu3 06+j/. 

In order to disentangle structural and Coulomb disor- 
der we compare isovalent (Eu,Nd) and heterovalent (Ca) 
substitutions. Based on the ionic radii^ 8 - - Y 3+ (101.9 
pm), Eu 3 + (106.6 pm), Nd 3+ (110.9 pm), Ca 2+ (112 pm) 
- we evaluate the cation mismatch via the radius stan- 
dard deviation, a r = 1.3,2.4,2.2 pm for Eu8%, Nd7.5% 
and Y .95Cao.o5Ba2Cu306+ y (Ca5%), respectively. Our 
data show that the Neel and superconducting transi- 
tions of the Nd and Eu substituted samples are identical 
within experimental error to those of the pure compound 
(Fig. [2]), whereas Ca5%, with the same oy as Nd7.5%, 
displays a markedly different behaviori 20 ' 41 We conclude 
that, for dilute Y replacements, isovalent disorder is by 
far less effective than heterovalent disorder. 

As a consequence the three families (Y100%, Eu8% 
and Nd7.5%) equivalently identify the cuprate system 
closest to the clean limit and in the discussion of the 
following subsections we shall refer generically to them 
as the YBCO system. 



IV. DISCUSSION 

A. Isovalent Y substitutions 

Disorder must be considered as an additional parame- 
ter controlling the phase diagram of oxides, distinguish- 
ing a clean and a dirty (heavily disordered) limit . 17 ' 18 The 
latter case is clearly represente d 14 ' 20 ' 35 by Yi_ x Ca 2 .Ba2- 
CuaOg+j, and by I^-zSr^CuC^. These examples con- 
sist of heterovalent cationic substitution in layers nearest 
neighbor to Cu02, determining a random distribution of 
Coulomb impurity potentials accompanied also by local 
structural distortions. The strong influence of the latter 
on the electronic properties^ 5 - is specifically demonstrated 
by the emergence of static stripes, evoked in La2-^Sr K - 
CUO4 by the La-Nd,E u 13 ' 46 partial substitution around 



B. Nature of the two magnetically ordered regimes 

Our data establish a coherent and systematic picture 
of the AF YBCO phase which is characterized by two 
distinct regimes, re-entrant and thermally activated, in 
agreement with previous indications 3 - 1 - on the magnetic 
behavior of La2-xSr x Cu04. The low temperature re- 
entrant regime witnesses the recovery of the full 2DHAF 
moment, fJ.2D, as it is shown in Fig [3^ (solid curve). This 
finding is reconciled with the opposite sign anomaly in 
the neutron magnetic Bragg peaks 7,8 considering that 
the neutron elastic intensity increases 8 - at low temper- 
atures along the 2D ridges (for q |j c*). Although the 
total elastic scattering does not allow a unique identifi- 
cation of the low temperature magnetic structure, the in- 
commensurate peaks 10 observed in La2-xSr x Cu04 above 
h c = 0.02 are tentatively interpreted in terms of a spiral 
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FIG. 6: (color on-line) Magnetic moments rap and tua 
(same as Fig. [3j compared with theory^ (dotted line) and 
experiments*- on La2-2,Sr x Cui_ z Znz04 (crosses). 
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FIG. 7: (color on-line) Scaling of the moment in the activated 
regime, itia, with Tn (for the solid line see the text). 



spin stated In YBa2Cu306+ y the experimental evidence 
of incommensurate order is missing, but similar mag- 
netic correlations are detected^ Stripes or spiral spin 
structures are indeed a viable candidates since they are 
locally collinear, in-plane and staggered, as required by 
the muon local field (Appendix [XJ and they can provide 
the large local moments that we detect. 

Another important feature of the re-entrant state is the 
narrow width of the observed moment distribution Amp, 
that reaches the value 0.4 only close to the critical point 
h c p, as shown in Fig.[5p. Notice that narrow widths and 
staggered in-plane collinearity extend also to the CSG 
phase coexisting with superconductivity, for h c < h < 
h c p. This is not what is expected of a canonical spin 
glass, where muons typically™ measure Am w 1 and 
randomly oriented local fields. 

The main signature of the CSG phase is the change 
from a rather large slope of the transition temperature, 
dT]\[/dh, below h c , to a very modest slope dTf /dh above. 
Notice that the Ta, the activation temperature defined 
by Eq. [2 shares this same modest slope, and the two 
temperatures, Ta and Tf follow the same curve in Fig. [5] 
This may be simply understood in terms of the disap- 
pearance of the high temperature, activated magnetic 
regime, above the critical density h c , while the mag- 
netic ground state detected in the re-entrant regime sur- 
vives as shown by the continuous behavior of mp{h) and 
ArriF{h) across h c , Fig. [31 This leads us to identify the 
same ground state in the CSG phase coexisting with su- 
perconductivity [h c < h < h c p) and in the re-entrant 
regime of the antiferromagnet, a common state that we 
propose to dub frozen antiferromagnet (FAF), in contrast 
with the thermally activated antiferromagnet, henceforth 
TAAF, established only above the crossover temperature, 
for h < h c . 

Insight on the nature of the FAF state is provided 
by Fig. [6] which shows that mp is actually following 
the initial slope typical of magnetic site dilution in 
cuprates^S (dotted line). A similar conclusion can be 
derived from recent low temperature muon results^ in 



G& x \ja.i.25Q&i.75-xC~u?,OG+y One way to look at this is 
to suppose that in the re-entrant regime holes are local- 
ized and each of them in first approximation cancels one 
Cu spin, leading to magnetic dilution with z = h. 

The dilution theory has recently founds excel- 
lent agreement with available La2- a; Sr a ;Cui-2Zn z 04 
dat a 51 ' 53 ' 54 (crosses in Fig. [5]), by taking into considera- 
tion a small additional frustrating effect of the local Zn 
impurity that slightly reduces the effective coupling. Also 
localized holes give rise to a frustration^ in this case of 
dipolar nature, that is often indicated as the cause of the 
dramatic reduction of Tjy(h). Our mp(h) data in Fig. |6l 
however, follow very closely the dotted curve of Ref. [50 
until the magnetic component disappears, showing that 
the leading phenomenon at low temperature is moment 
dilution and not frustration. This result indicates that 
holes, either by incoherent freezing or by self-organization 
produce only a rather moderate frustration in the AF 
background of the FAF state. 

Conversely the slope of the thermally activated mo- 
ment vriA{h) is much steeper than the predictions of the 
dilution theory, even at extremely low doping, as the 
dashed and dotted curves show in Fig. O It is noteworthy 
that the vanishing of tua, the divergence of Am a (Fig. [3]) 
and the suppression of TV (Fig. [2J , all occur at the same 
critical density h Cl thus linking all these properties to the 
specific TAAF phase alone. 

Figure [JJ shows that the activated moment rriA(h) 
and Tjv(/i) scale roughly linearly with each other over a 
large range of h (the solid line is obtained by plotting 
the previously denned functions m^/i, oo) vs. T^{h) 
with h as an implicit parameter). It implies that both 
niA and Tn scale with an effective exchange coupling 
J e //(/i), strongly reduced by the activated process, either 
by finite-size effects 3 ** or by a very large frustration,—^ 
much larger than that present in the re-entrant regime. 

The magnetic moment in the TAAF state follows the 
phenomenological fit of Eq. [5] Nominally, this fit justifies 
the low temperature regime in terms of hole localization. 
Here the muon is directly sensitive only to spin dynam- 
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ics, hence, recalling that spin/charge separation may take 
place in cuprates^ we cannot actually prove directly that 
charges freeze as well. Freezing processes with similar low 
energy scales are detected also by /xS R 14 ' 34 and NQ R 36 i 58 
T^ 1 relaxations (the latter being sensitive also to charge 
dynamics) in La2- x Sr x Cu04 and in YBa2Cu306+ y as 
well. 

Direct evidence that thermal activation regards also 
charge degrees of freedom comes from transport, where 
low temperature variable range hoppin g 15 ' 59 is detected 
both in dirty La2- :E Sr ;r Cu04 and in clean YBa2Cu306+ y , 
even at very low hole densities^ The energy scale of 
charge activation is comparable to the one we observe 
here (resistivity at the lowest doping in Ref. 38, = 380 
K, shows a low temperature divergence starting below 70 
K). The TAAF phase is therefore the same that strongly 
violates the Mott-Ioffe-Regel limit (Fermi wave vector 
times mean free path much less than one), but still dis- 
plays a metallic transport character— even at very low 
hole densities. The coincidence of the onset of this so- 
called bad metaMl state, with the reduction of the mag- 
netic moment around Ta consolidates the association of 
thermal activation with charge carriers implied by Eq. [2] 
This point deserves a systematic investigation by trans- 
port in very low doped single YBa2Cu306+ y crystals, 
scarcely studied to date. 

Notice that when the TAAF phase disappears, for h > 
h Cl the samples become electronically inhomogeneous, 
since the well known nanoscopic phase separation 16 takes 
place. In this case the bad metal state exists already at 
T = in the percolating superconducting background, 
but inside the interspersed nanoscopic magnetic clusters 
it develops only above Tf. 

If thermally activated holes determine directly the re- 
duction of m in the TAAF state they must be responsible 
of a much larger frustration than independently localized 
holes. Such a view could be reconciled with recent ex- 
tensions of the spiral model 9 -, where increasing hole local- 
ization lengths, are shown to reduce the staggered 
moment very effectively. The same approach successfully 
accounts for activated hole transpor t 61 ' 62 , i.e. for the ef- 
fect of spin texture on carrier mobility, but none has yet 
taken in due consideration the reverse effect of hole mo- 
tion on spin texture, which is very relevant, as we show. 
The two-regime-behavior of Fig. U would then suggest 
that the localization length is not temperature indepen- 
dent in cuprates: while a rather large length k^ 1 (cfr. the 
value of 100 pm quoted in Ref [9j for La2- x Sr x Cu04) is 
appropriate for the thermally activated regime, a smaller 
length would be characteristic of the FAF state. Such 
a temperature-dependent localization scale suggests that 
at very low temperature the effects of charge localiza- 
tion in a lightly doped Mott-Hubbard insulator are not 
fully captured by the standard impurity model of doped 
semiconductors investigated so far in the literature. 

However, as we have noted below Eq. [SJ the fit of the 
muon measurements to this function by itself does not 
rule out alternative stripe based ideas. For instance a 



similar connection with thermally activated hole density 
could emerge considering stripe domains of finite length 
and with rather large d — 1/2 ft, separation^ 3 - at the low 
densities considered here. In these conditions metallic ID 
stripes would not percolate across the domains, whereas 
the antiferromagnetic regions would. Thermally acti- 
vated hopping would then switch the system to the op- 
posite condition of percolation for transport and separa- 
tion of independent magnetic clusters, at the same time. 
This process is already favored by the reduced effective 
exchange across the metallic ID stripe, and it would lead 
to finite-size reduction of both m and TV- 

Summarizing, both spirals and stripes could qualita- 
tively justify the YBCO /xSR results, but only further 
theoretical work can determine if either of these hypothe- 
ses reproduces the observed crossover and the associated 
large change in the static local moment. 

Another important result of the present work is that 
the critical value, h c — 0.056(2), found for the suppres- 
sion of Tjv(ft) and rriA{h) coincides with the critical value 
for the onset of T c (h), h s — 0.056(2). Both values agree 
within our errorbars with a recent very careful single 
crystal determinations^ of h s — 0.0539(1). According 
to our data also Tjv(/i), in the YBCO clean limit, fol- 
lows a parabolic behavior analogous to tha t 65 ' 66 of T c {h) 
and very different from theoretical predictions. 56 Notice 
that older systematic results, 5 " 7 reported only vs. oxygen 
content, without hole doping calibration, could not pro- 
vide such information. Therefore our results indicate 
that if superconductivity and the TAAF state were the 
only competing orders, h c — h s would represent a true 
quantum critical point (QCP). The critical density is 
rather lower than that more often indicated for a QCP in 
cuprates^ 7 -, but other authors^ have considered a similar 
additional location. The real-world system however does 
not develop the QCP, because it switches instead to the 
FAF ground state throughout a wide doping range, from 
h fa to well past the onset of superconductivity. Ad- 
ditional evidence on this point comes from our parallel 
work on the heterovalent substitutions^ 



V. CONCLUSIONS 

We have extensively investigated YBCO in the anti- 
ferromagnetic low doping region as a model for the clean 
limit cuprate. Since the influence of disorder induced by 
Nd and Eu is found to be negligible, we extract the pa- 
rameters of the low temperature re-entrant and of the 
thermally activated regimes from the whole set of data. 
The crossover between these two regimes is undoubtedly 
associated with the closely related thermal activation of 
both spin and charge degrees of freedom. This is imple- 
mented in a single phenomenological model for the tem- 
perature and doping dependence of the moment m(h, T), 
which fits our data throughout the entire investigated 
range. 

Our analysis leads to the identification of a common 
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ground state for very different doping regimes at the two 
sides of the metal-insulator transition: the so-called Clus- 
ter Spin Glass phase and the re-entrant antiferromagnetic 
phase. We dub this unique state a frozen antiferromag- 
nct. 

We have shown that in the frozen antiferromagnetic 
state the ordered magnetic moment is subject to negligi- 
ble frustration and it follows the same trend of magnetic 
site dilution. Conversely the fast drop of T^{h) and of 
the staggered moment, m^(/i) are characteristic of the 
thermally activated antiferromagnetic phase. 

Finally, our data indicate that the main boundaries in 
the clean-limit phase diagram, describing the two main 
order parameters, magnetic and superconducting, both 
follow parabolic curve o 16 i 41 i 66 ' 69 sharing a zero temper- 
ature intersection at h = h c = h s . At face value this 
could mean that h c = 0.056 is a quantum critical point 
for the cuprate clean limit. The QCP is then superseded 
by the FAF behavior of the real-world compound at low 
temperatures. 



Appendix A: The muon local field 

Three types of muon stopping sites are distinguished in 
cuprates. The dominant one, observed in all compounds, 
irrespective of both the details of the perovksite structure 
and doping, is bound to the apical oxygen (AO), just 
above or below the Cu02 plaquette (inset of Fig (5]). The 
low temperature value of the field probed by muons at 
this site is B L AO = 40 mT for La 2 Cu0 4 , Ref. [zS and 
B Y AO = 30 mT for YBa 2 Cu 3 6 , Ref. i, respectively. A 
second site is specific of YE^CusOg+j, for y > 0.2 and 
it is attributed to muons bound to chain oxygen (CO), 
with roughly half the field intensity. A third site, directly 
bound to the oxygen ions of the plaquette (PO), with 
large internal fields, is observed only£i in undoped Y- 
Ba 2 Cu 3 6 . 

The AO site is the main subject of our analysis, and its 
assignment is based on numerous observations j22r— One 
of them is that the value of B \ Q at zero temperature and 
doping is reproduced by simple dipolar sums: if S(h,T) 
is the average staggered Cu spin defined in Sec. IIII1 Si is 
the spin direction at site i, Ti = is the vector joining 
the muon site to the i-th Cu spin and g is the Lande 
factor, we can write 



B A o(h,T) 



47T 



E 



3(Sj 



\D\S(h,T), 



S(h,T) 
(Al) 



The zero temperature value is reproduced in this site 
with the standard — 2D spin reduction, gS(0, 0) = 0.6. 
A small isotropic super-hyperfine contribution cannot be 
a-priori excluded, but it would preserve the validity of 
Eq. lAll with the simple addition of a term AS\ to D, Si 
being the Cu spin nearest to the muon. 
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FIG. 8: (color on-line) YBa 2 Cu306+ H , y = 0.27, h = 0.023: 
FFT amplitude of the precessing asymmetry a three tempera- 
tures, vertically shifted for clarity. Inset: muon AO site, with 
apical oxygen and the CuQ2 plaquette. 



The sums defining vector D in Eq. lAll converge rapidly, 
thanks to the alternating sign of Si in the AF structure 
(the staggered moment). Even in the very underdoped 
superconductor YBa2Cu3 06.35, where in-plane magnetic 
correlation lengths are shortest^ they never fall below 
£a t b 12 A, and the relative contributions to Eq. I All for 
fj > £a,&/2 are only a few percent. Hence the local field 
may be written as a constant, D = \D\, times S(h,T) 

The rapid convergence of the dipolar sums and the 
tensorial nature of the interaction guarantee that D does 
not change as long as the local spin directions Si, nearest 
neighbors to the muon, remain i) collinear, ii) parallel 
to the ab plane, and ii) staggered. These points follow 
from the consideration that the value of Bao obtained 
in the simple Neel arrangement is actually invariant un- 
der in-plane Si rotations for a muon site coinciding with 
that of the apical oxygen. The AO site is displaced only 
by 110 pm from this symmetry position and the invari- 
ance is roughly preserved. Conversely if either the mo- 
ment were drastically reduced, or the local spin arrange- 
ment did change significantly from staggered, collinear 
and in-plane, the low temperature value of the local muon 
field of Eq. IA1I would also change. For these reasons we 
conclude that the nearly constant local field BAo{h,0) 
(rescaled to mp{h) in Fig. [3Ji) requires both nearly con- 
stant D and S. An unvarying D in turns implies that a 
collinear staggered in-plane spin arrangement is at least 
locally preserved at low temperature for all hole densities. 

The argument discussed above could in principle be 
questioned in the presence of an incommensurate mag- 
netic structure with long, /i-dependent wavelength, giv- 
ing rise to broad distributions of local fields with a sharp 
peak at the Van Hove singularity. Since only this peak 
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FIG. 9: (color on-line) Lattice parameter c from XRD versus 
hole density from thermopower. The solid line is the calibra- 
tion from Ref. l66l . 

would be observed and fitted, any doping-dependent 
change could modulate not only the peak position, but 
also its spectral weight (i.e. reduce the fitted initial union 
asymmetry). However we always detect an unchanging 
spectral position and the full initial muon asymmetry, 
ruling out this possibility. 

Let us conclude by reviewing the experimental evi- 
dence: the negligible variation of the low temperature 
AO field throughout the magnetic phase diagram is am- 
ply discussed in Sec. IIIII The field distribution itself 
does not change drastically as it is witnessed by the very 
similar low temperature asymmetries in the six panels 
of Fig. HH spanning from h = 0.02 (top) to h = 0.07 
(bottom) . If local correlations did vary significantly we 
would observe the appearance of distinct frequency com- 
ponents, signalling different internal fields, whereas the 
only notable change is that the damping of the preces- 
sions (i.e. the width of the field distribution) increases 
with hole content, as shown also in Fig. \3]p. 

The field distribution does not change appreciably with 
temperature either. Figure [8] shows representative Fast 
Fourier Transform (FFT) spectra of the muon precessing 
asymmetry below, around and above Ta- They display 
two peaks corresponding to the AO and CO sites, which 
shift to lower frequency as S(T) is reduced, preserving 
their shape. 



Appendix B: Experimental details 

High quality polycrystalline samples of YBa2Cu3 06+i,, 
Yo.925Ndo.o75Ba 2 Cu30 6+ j / and Y .92Euo.o8Ba2Cu 3 6+ , / 
with 0.05(2) < y < 0.48(1) are obtained by sintering. 
The YI00% series is the same of Ref. The absence of 
spurious phases is confirmed by systematic x-ray Rietveld 
refinements, and the oxygen content has been fixed by a 
topotactic technique, which consists in the oxygen equi- 
libration of stoichiometric quantities of two end mem- 
bers (y = 0.05(2) and y = 0.98(1)), packed in sealed 
vessels with minimum free volumei 76 ' 77 Low tempera- 



ture annealing produces homogeneous samples, with well 
ordered CuO chains in the basal plane, with an abso- 
lute error of Sy — ±0.01 in oxygen content per formula 
unit (cross checked by iodometric titration and thermo- 
gravimetry on each sample). The mobile carrier con- 
tent is obtained from systematic measurements of the 
room temperature Seebeck coefficient, 78 with a calibra- 
tion procedure^ of the exponential dependence of S vs 
h, performed on empty-chain Yi^Ca^I^CusOe sam- 
ples. Figure [9] shows that this calibration agrees very 
well with that extracted^ from the lattice parameter c. 
Notice that our determination of a common hole doping 
density scale for all compounds is not based on the a- 
priori assumption of a common equal doping density for 
optimal superconductivity h Q rs 0.16, as it was done for 
instance also in our previous worki 16 ' 79 . Slight deviations 
between Fig.[2]and these previous published Y100% data 
are due to the different calibration. 

Superconducting behavior characterizes samples with 
h > h s = 0.056, whose critical temperature T c is deter- 
mined from the linear extrapolation of the 90% to 10% 
drop in the diamagnetic susceptibility measured in a field 
of fi H = 0.2mT. 

The /iSR experiments were performed on the MUSR 
and EMU spectrometers at the ISIS Facility of the 
Rutherford Appleton Laboratory (UK) in longitudinal 
geometry with zero external field. 80 



Appendix C: Zero-Field /uSR analysis 

The positron emission in the muon decay is asymmet- 
ric, since it is preferentially directed along the muon spin 
I. For a polarized muon ensemble the measured asym- 
metry is thus directly proportional 80 to the average spin 
self-correlation function, A(t) = A I z (0)I z (t), where z is 
the axis of the detector and Aq an experimental param- 
eter subject to calibration. 

Typical asymmetry functions in magnetically ordered 
and paramagnetic phases are shown in Fig. [TO] The left 
column refers to Y100%, the right one to Nd7.5% (very 
similar data are obtained on Eu8% as well). From top 
to bottom three different samples for each class of com- 
pounds span the magnetic phase diagram, from very low 
to intermediate doping (h sa 0.02, h « 0.04, and h « 0.07 
respectively). In each panel the time evolution of the 
muon asymmetry is shown for three temperatures, the 
lowest, one close to the magnetic transition and one just 
above it. 

The self-correlation function yields a simple relaxing 
term, if the internal magnetic field B^ is either vanish- 
ing or longitudinal to the initial polarization direction 
1(0), and an additional precessing term, if B^ is trans- 
verse. Therefore in all six panels the best fit of the sam- 
ple asymmetry at the lowest temperature, well below the 
magnetic transition, show both a simple relaxation and 
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FIG. 10: (color on-line) Muon asymmetry above Tjv, just be- 
low and for T <C Tjv. Left: three YBa2Cu30e+i/ samples, 
right: three corresponding Yo.92sNdo.o75Ba2Cu306+; / sam- 
ples, from nearly undoped (top) to above the onset density 
for superconductivity. 



precessions, following 



A(t) = A 



f, 



cos + / 



i-\t 



, (CI) 



where 7 = 0.8514 • 10 9 (sT)" 1 is the muon magneto- 
gyric ratio, /*' are the transverse and longitudinal muon 
fractions, experiencing relaxation rates di and A, respec- 
tively. The index i spans the two muon sites, AO and 
CO (Appendix [AJ , where muons probe two distinct lo- 
cal fields Bao an d Bco ■ The longitudinal term does not 
distinguish the two sites. The Gaussian form of the trans- 
verse relaxation indicates a distribution of static internal 
fields p(Bi), whose second moment (AB^ 2 ) is connected 

to the rate by <x; = 7 (AB,; 2 ) 1 ' 2 . The experimental data 
generally include also a small non-relaxing term due to 



muons stopping outside the sample. 

A simple geometric argument shows^S that for a sin- 
gle site in a magnetically ordered polycrystal the experi- 
mental muon fractions correspond to an average over all 
crystal orientations that yields / = 1/3 and /* = 2/3 
(in our case /' must be replaced by f AO + fco)- When 
magnetic and non-magnetic phases are present in the 
sample, since muons are evenly implanted in both, their 
fractions reflect the relative volumes of the two phases. 
The relative magnetic volume v m = V m /Vtot is given by 
V m = 3(1 — / )/2 = 3/*/2, which is the quantity shown in 
Fig. [1}d. Its smooth drop across the magnetic transition 
temperature signals some oxygen inhomogeneity, produc- 
ing a Gaussian distributio n 16 ' 34 of mean T m and width 
AT m , hence the relative magnetic volume can be fitted 
to the function v m {T) = [1 - erf((T - T m )/ v / 2AT m )]/2. 

The middle plot in each panel of Fig. [TO] shows data 
just below T m characterized by over-damped precessions 
and increased value of the longitudinal fraction / . Best 
fits allow on ly for one undistinguished transverse fraction, 
with B = VAB 2 . 

At the highest temperatures, well above the magnetic 
transition, samples are fully in the paramagnetic state 
and the sample asymmetry is fitted to the function 



A(t) = A e 



(C2) 



where a n is the static contribution to the relaxation rate 
arising from neighboring nuclear moments. 

Above T « 250 K muon diffusion sets in, preventing 
the direct measure of the transverse precessions. For this 
reason high temperature points are missing in Figs. [TJt 
and m However the Neel temperature can still be ob- 
tained from the longitudinal fractions, 34 which are not af- 
fected by the incipient muon diffusion, as long as an over- 
damped transverse fraction can still be distinguished. 
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